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1. Introduction 

The last few years have seen remarkable progress in the numerical study of QCD at finite 
chemical potential }X. Still, the various methods [|l|, ^ ^, ^ suffer from systematic uncertainties, 
which limit their range of reliability to about 7 ^ 1.0. For a recent review, see Ref. [^. We try to 
address this apparent limitation by using a canonical approach 0], where we focus on the matter 
density p, rather than the chemical potential. The method is particularly appropriate to explore few- 
nucleon systems at low temperature, and in principle, allows to study the bulk properties of nuclear 
matter and the nuclear interactions. Here, we extend its use and determine the phase boundary 
between the confined phase and the quark gluon plasma, as illustrated in Fig.|l|. 




^^ p 



Figure 1: Sketch of the conjectured QCD phase diagram in the grand-canonical and canonical formalism. 

The left figure shows a sketch of the conjectured phase diagram of QCD in the T-^ plane. 
At small chemical potential, the phase transition is a rapid crossover, which ends in a second 
order endpoint, followed by a first order transition. Correspondingly, the right figure illustrates the 
transition in the T - p plane. The first order transition is manifested in a co-existence region. In 
this proceedings, we describe how we identify the co-existence region and how we determine the 
phase diagram in the T-pL as well as in the T-p plane. 

2. Partition Functions 

We construct the canonical partition function Zc{T,Q) by fixing the number of quarks N = 
J cPx \j/{x) 7o to Q- We insert a 5-function in the grand canonical partition function Zgc{T,IJ.) 

ZciT,Q) = I [DU][D^][DW] e-sAU-T]-sAu.^.'i'J,f,]g . (2.1) 

The 5-function admits a Fourier representation 8 (^N — Q) = J dp.; e'^'^'^^S) y^/e recognise //i/ = 
ijljT as an imaginary chemical potential and exploit the ^^-periodicity [^] in jLt/ of Zgc{T, }X = ifXi) 

Zc{T,Q) = ^ r dfiie-'Qf''ZGc{T,ifijT) = ^ T d (^) e-''""^ ZcdT J^i,) . {12) 

Thus, the canonical partition functions are the coefficients of the Fourier expansion in imaginary /i 
of the grand canonical partition function. As a consequence of the ^^-periodicity, the canonical 
partition functions are zero for non-integer baryon number B = Q/3. 

From the canonical partition functions, the grand canonical partition function can be recon- 
structed using the fugacity expansion (in fact a Laplace transformation) 

Zgc{T,p) = r dpe'''PrZciT,p)= r dpe-r(n^'py'^P^ (2.3) 
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with the baryon density p = f and the Helmholtz free energy density f{T,p) = —ylogZc{T,p). 
The relation between baryon density and chemical potential can be expressed as (p)(/x) or ^ii (p): 



(P)(M) 



1 



■Jdppe'''PrZc{T,p) 



or 



m(p) 



1 dfip) 
3 dp 



(2.4) 



While the first expression is exact in any volume, the second is obtained via a saddle point approx- 
imation (exact in the thermodynamic limit) and may have more than one solution when solving for 
the baryon density at a given chemical potential, see Fig.^. We discuss this issue in detail below. 

3. Method 

Following [^], we express the canonical partition function in a ratio, which can be measured 
by Monte Carlo simulation as an expectation value: 



Zc(B,j8) 



1 



■ / [DU] . 



^2n 



ZGc(i3o,'M/o) 

-iSB^ det([/;//i/) 



det{U;iiiio 



1 



det(?7; //!/(, ) 



)j3c 



In J-iz 
Zc{U-B 



det(?7;//i/n 



det(?7;//i/) 
det{U;i}Xi^) 

(3.1) 



where Zcci^o^ip^h) is the grand canonical partition function sampled by ordinary Monte Carlo 
methods, here for notational simplicity at jSq = /3. The Zc(?/;S)'s are the Fourier coefficients of 



the fermion determinant for a given configuration {U}. Although the average in Eq.(|3.1|) should be 
real positive, the individual measurements are complex, with a sometimes negative real part. This 
is how the sign problem manifests itself in our approach. Moreover, a reliable estimate depends 
on a good overlap of our Monte Carlo ensemble with the canonical sector B at temperature j8. We 
address this issue by following the idea of Ref. [|I|] and including both confined and deconfined 
configurations in our ensemble. Indeed, we supplement the ensemble at (j3f(/i = 0),/^ =0) with 
additional critical ensembles at imaginary chemical potential, non-zero isospin chemical poten- 
tial, and ensembles generated with an asymmetric Dirac coupling [Q] - in principle, any ensemble 
is allowed. We then combine all this information about a particular canonical partition function 
Zc(B,/3) by Ferrenberg-Swendsen reweighting JToj]. 

The Fourier-coefficients of the determinant Zc{U ; Q) are calculated exactly [Q]. In the tempo- 
ral gauge (?74(x, t) = \ except for t =Nt — I), the staggered fermion matrix M in the presence of a 
chemical potential can be written in the form 



Bo 
1 





1 ... 

Bi 1 
-1 B2 1 



ul_ 



-llaN,\ 



\^-f/;v,-ie^''^' ... -1 %_i ) 




Un,-i 



where the B,'s contain all space-like contributions. Ref. [ |1 IQ showed that the determinant can be 
computed for any chemical potential at the cost of diagonalising the so-called "reduced matrix" P. 
The determinant is given in terms of P's eigenvalues Ai, . . . , Agy, where V is the spatial volume: 



6V Q=3V 

detM([/; H)=e^^ H + e'^"^') = L Zc{U;Q)e-Q^"^' 

i=l Q=-3V 



(3.2) 
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Matching term by term, we then solve for the Fourier coefficients Zc{U;Q). This delicate step 
requires a special multi-precision library. The diagonalisation of the reduced matrix P is computa- 
tionally intensive and takes ff{V^) operations. 



4. Results 

We study the Helmholtz free energy F{B) = — Tlog^j^ in a theory of four degenerate 
flavours of staggered quarks with mass ma = 0.05 {'f = 0.2, nij^ 350 MeV) on a small 6^ x 4 
lattice with volume ~ (1.8fm)^. For the quark mass we chose, the phase transition is first order at 
jJ. =0, and presumably remains first order for all chemical potentials. 

We "scan" the phase diagram by varying the baryon density at fixed temperature, see Fig.^. 

We measure ^^j. and assume the validity of the saddle point approximation to equate this 

quantity with following Eq. (2.4). This assumption will be tested in Fig.^. 
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Figure 2: (left) A sketch of the "scans" in our phase diagram, (right) The derivative of the free energy at 
fixed temperature as a function of the baryon number (baryon density). In the saddle point approximation, 
the y-axis is ^ = ^ ' -. 

Note that accurate results are obtained, up to high densities (> 5 baryons/fm^) and large 
chemical potentials (^ ~ 2). The first order phase transition and the associated metastabilities 
are clearly visible in the "S-shape" of jiB). The low-density regime can be reasonably well de- 
scribed by a simple hadron resonance gas Ansatz, jj = 3f{T) sinh(3^) with f{T) as the only free 
parameter. The high-density regime almost corresponds to a gas of free massless quarks jj = 



Nf (^) + (^j) when taking cut-off corrections [12] into account. The solid line in Fig.||(right) 
is obtained by fitting the linear and cubic terms in this expression. Instead of the free value 1, the 
fitted coefficients are 0.82(2) and 1.94(6) respectively. Thus, the equation of state for the quarks in 



the plasma phase differs little from the Stefan-Boltzmann law. This has been observed also in [ 13] 



and JlH]. Note that we find the same jjij) dependence in the plasma phase at all temperatures. 

For a given temperature T, we identify the boundaries pi and p2 of the co-existence region 
and the critical chemical potential /i as follows. Equality of the free energy densities in the two 
phases, /(pi) - 3;Upi = /(P2) - 3;Up2, implies 



r'dp{f'{p)-3^)=0. 



(4.1) 



Since /'(p) is the quantity measured in Fig.^, we determine Pi,P2 and /i by a "Maxwell construc- 
tion" illustrated in Fig.^ (left) for the temperature y = 0.92. The value of j defining the horizontal 
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line is adjusted to make the areas of the two "bumps" in the S-shape equal.' The two outermost 
crossing points define p\ and p2, the boun( 
the value of the critical chemical potential. 



crossing points define p\ and p2, the boundaries of the co-existence region. Here, ^ = 1.06(2) is 
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Figure 3: (left) The Maxwell construction allows to extract the critical chemical potential and the bound- 
aries of the co-existence region, (right) Comparing the saddle point approximation (red) with the fugacity 
expansion (blue). Strong finite-size effects in the latter obscure the first-order transition. 



We can cross-check this result by making use of the fugacity expansion Eq.(2^), see Fig.|| 
(right). For a given chemical potential, we measure the baryon number (B)(/i). We see a jump at 
the same value ^ w 1 .06, but the rounding due to finite size effects is very strong. In contrast, our 
criterion for criticality (equality of the free energies) has exponentially small volume corrections. 

In Fig.0 we present the phase diagrams in the T-jx as well as in the T-p plane. On the left, we 
summarise results from various methods, all for the same theory: 4 flavours of staggered quarks 
with am = 0.05, Nt = 4 time-slices; only the spatial volume varies as indicated. We have repeated 
(blue) the study of [|l|] (green), using multi-parameter reweighting on one ensemble generated at 
(j8c,/x = 0). We identify the phase transition via the peak of the specific heat instead of Lee- 
Yang zeroes, and obtain consistent results. However, the "sign problem" dramatically grows with 
increasing chemical potential, as shown by the average sign in the figure. Moreover, our statistical 
en^or, based on jackknife bins as in [|l]], does not reflect the true inaccuracy. 

The parabolic fit [^] is consistent with the black points [Q]. Both methods perform an analytic 
continuation from imaginary jj., for which the systematic errors are hard to quantify. Our new 
results are shown in red. There is no strong inconsistency with other results, but we observe a clear 
sign of bending down starting at ^ ~ 1.3. In fact this must happen, if the critical line is to reach 
the value ajXc = 0.35 at /3 = 0, predicted from a strong coupling analysis [15]. In the T-p plane, 
the densities at the boundaries of the co-existence region seem to remain constant for T < 0.857^ 
already, with Pqgp = 1.8(3)B/fm^ and Pconfined = 0.50(5)B/fm^. The latter is a plausible value for 
the nuclear density in our 4-flavour, = 350 MeV QCD theory. 



5. Conclusions 

We study QCD in a canonical framework, which is promising for the study of few-nucleon 
systems at low temperature, but proves also capable of exploring high density regimes {pL/T <2) 

'The area of each bump gives the free energy required to build two planar interfaces. The corresponding interface 
tension is ^ ~ 35 - 45 MeV. 
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Figure 4: (left) The phase diagram in the r-/i -plane, (right) The phase diagram in the T-p -plane. 



at temperatures T > 0.87^. . We have determined the phase boundary between the confined phase 
and the quark gluon plasma in both the T-p and the T-jX plane. In the latter, our results are in 
agreement with the literature, however we observe a bending down of the critical line at ^ ~ 1.3. 
The two phases can be rather well described by the hadron resonance gas at low densities and by a 
weakly interacting massless gas at high densities. 

We thank the Minnesota Supercomputer Institute for computing resources. 
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